F/G  20/3 


AD-A107  832 

UNCLASSIFIED 


MOORE  SCHOOL  OF  ELECTRICAL  ENGINEERING  PHILADELPHIA  PA 
HIGH  RESOLUTION  3-D  WAVELENGTH  DIVERSITY  IMAGING. (U) 

SEP  81  N  H  FARHAT  AFOSR-BO-0281 

E0/M0-4  AF0SR-TR-B1-0731  NL 


UNIVERSITY  of  PENNSYLVANIA 
The  Moore  School  of  Electrical  Engineering 

PHILADELPHIA,  PENNSYLVANIA  19104 


Approved  for  public  release 

distribution  unlimited. 


gj  11  10  0o8 


UNIVERSITY  OF  PENNSYLVANIA 
THE  MOORE  SCHOOL  OF  ELECTRICAL  ENGINEERING 
ELECTRO-OPTICS  ANN  MICROWAVE-OPTICS  LABORATORY 
PHILADELPHIA,  PENNSYLVANIA  19104 


FINAL  REPORT 

HIGH  RESOLUTION  3-D  WAVELENGTH  DIVERSITY  IMAGING 


Prepared  for  the 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH  N/E 
BUILDING  410,  BOLLING  AIR  FORCE  BASE 
WASHINGTON,  D.C.  20332 


Under  Grant 
AFOSR-6 0-0281 
Covering  the  Period: 
8/1/80  to  7/31/81 


by 

N.  H.  Farhat 


September  25,  1981 


EO/MO  Report  No.  4 


r  - 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Detm  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
!  BEFORE  COMPLETING  FORM 

|1.  REPORT  NUMBER 

AFOSH-TR.  81  -073  1 

m 

3.  RECIPIENT'S  CATALOG  NUMBER 

J  5-23905 

4,  TITLE  (end  Subtitle) 

HIGH  RESOLUTION  3-D  WAVELENGTH  DIVERSITY 

IMAGING 

5  TYPE  OF  REPORT  6  PERIOO  COVERED 

Final 

8/1/80  to  7/31/81 

8.  PERFORMING  O^G.  REPORT  NUMBER 

EO/MO-4 

7.  author^; 

N.H.  FARHAT 

6.  CONTRACT  OR  GRANT  NUMBERfvJ 

AF0SR-8 0-0281 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

University  of  Pennsylvania 

The  Moore  School  of  Electrical  Engineering 

200  S.  33rd  St.,  Phila.,  Pa.  19104 

to.  PROGRAM  ELEMENT.  PROJECT,  TASK 

AREA  6  WORK  UNIT  NUMBERS 

6  // *-2- 

M.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

United  States  Air  Force 

Air  Force  Office  of  Scientific  Research 

Building  410,  Bolling  A.F.B.,  D.C.,  20332 

t2.  REPORT  date 

September  25,  1981 

1 3.  NUMBER  OF  PAGES 

56 

u.  MONITORING  AGENCY  NAME  4  ADDRESS!'//  different  from  Controlling  Office) 

15.  SECURITY  CLASS,  (of  thie  report) 

15*.  DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  ( of  thi»  Report) 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (oi  the  mbatrect  entered  In  Block  20,  It  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 

Three-dimensional  imaging,  wavelength  diversity,  polarization  diversity, 
tomography,  projection  hologram,  target  derived  reference,  aperture 
synthesis. 


19.  KEY  WORDS  (Continue  on  teveree  tide  if  neceeemry  end  Identify  by  block  number) 


20.  ABSTRACT  (Continue  on  reeeree  tide  If  neceeemry  end  identify  by  block  number) 

physical  optics,  vector  formulation  of  microwave  imaging  of  perfectly 
conducting  objects  by  wavelength  and  polarization  diversity  is  presented. 
The  results  provide  the  theoretical  basis  for  optimal  data  acquisition  and 
three-dimensional  tomographic  image  retrieval  procedures.  These  include: 
(a)  the  selection  of  highly  thinned  (sparse)  receiving  array  arrangements 
capable  of  collecting  large  amounts  of  information  about  remote  scattering 
objects  in  a  cost  effective  manner  and  (b)  techniques  for  3-D  tomographic 


^  Ui  (  X  ■  V  jfi  *  ^  . 

SECURITY  CLASSIFICATION  dF  THIS  PAGE  (When  Dmie  Entered) 


00  1  JAN ^73  1473  EOITION  OF  I  NOV  65  IS  OBSOLETE 


SCCUWTV  CUASSIFICATIQ^or  THIS  PAGEfWhm  Dmtm  Bnfrnd)  _ 

A 

Image  reconstruction  and  display  in  which  polarization  diversity  data 
is  fully  accounted  for.  Data  acquisition  employing  a  highly  attractive 
AMTDR  (Amplitude  Modulated  Target  Derived  Reference)  technique  is  dis¬ 
cussed  and  demonstrated  by  computer  simulation.  Equipment  configuration 
for  the  implementation  of  the  AMTDR  technique  is  also  given  together  with 
a  measurement  configuration  for  the  implementation  of  wavelength  diversity 
imaging  in  a  roof  experiment  aimed  at  imaging  a  passing  aircraft. 

Extension  of  the  theory  presented  to  3-D  tomographic  imaging  of 
passive  noise  emitting  objects  by  spectrally  selective  far  field  cross¬ 
correlation  measurements  is  also  given. 

Finally  several  refinements  made  in  our  anechoic -chamber  measurement 
system  are  shown  to  yield  drastic  improvement  in  performance  and  retrieved 
image  quality. 


i 


('  jL  >  <• — -JL 


SECURITY  classification  o 


AGEfWhwn  D*r»  Er> 


The  findings  in  this  report  are  those 
of  the  author  and  are  not  to  be  inter¬ 
preted  as  the  official  position  of  the 
Air  Force  Office  of  Scientific  Research 
or  the  U.S.  Government . 


TABLE  OF  CONTENTS 


Page 

1.  Introduction  1 

2.  Accomplishments  3 

2.1  Vector  Formulation  of  3-D 

Tomographic  Imaging  by  Wavelength 
and  Polarization  Diversity  4 

2.2  Target  Derived  Reference  Work  20 

2.3  Three  Dimensional  Imaging  of 

Incoherently  Radiating  Objects  29 

2.4  Hybrid  Data  Processing  and  3-D 

Image  Display  33 

2.5  Design  of  Roof  Experiment  37 

2.6  Improved  Measurement  System  and 

Experimental  Results  39 

3.  Conclusions  47 

4.  References  51 

5.  Publications  52 

6.  Personnel  53 


7. 


Appendix  I 


54 


1. 


INTRODUCTION 


The  motives  for  the  use  of  electro -magnetic  radiant  energy  in  the 
microwave  range  in  the  formation  of  images  of  scattering  or  radiating 
objects  are  numerous  and  varied.  Significant  among  them  are  the  ability 
of  microwaves  to  penetrate  through  inclement  weather,  smoke,  dust  and  a 
variety  of  solid  lossy  dielectrics  and  the  availability  of  a  sophisticated 
microwave  technology  that  permits  automated  and  accurate  direct  measure¬ 
ments  of  the  amplitude  and  phase  of  the  object  scattered  or  radiated  fields 
over  a  wide  range  of  frequencies.  A  conventional  approach  to  microwave 
imaging  based  on  scaling  of  camera-like  imaging  concepts  from  the  optical 
domain  to  the  microwave  range  is  impractical  because  of  the  cumbersome 
size  and  unacceptably  high  cost  of  the  microwave  apertures  required  to 
produce  high  resolution  images  especially  when  the  object  is  remote.  There 
is  also  the  problem  of  rapid  deterioration  of  longitudinal  resolution  with 
range  which  prevents  the  attainment  of  good  depth  resolution  in  the  image. 
Attempts  to  scale  optical  holography  concepts  to  the  microwave  range  suffer 
from  the  same  constraints  with  the  additional  problem  associated  with  wave¬ 
length  scaling  that  occurs  when  the  microwave  hologram  is  scaled  down  to 
the  optical  domain  in  order  to  view  the  reconstructed  image  in  laser  light 
as  in  an  optical  hologram.  The  wavelength  scaling  problem  prevents  the 
viewing  of  a  true  3-D  image  of  the  object  as  in  optical  holography  even 
when  the  object  is  not  so  remote  and  a  Fresnel  hologram  has  been  recorded. 
Furthermore  the  same  speckle  noise  which  plagues  most  coherent  optical- 
imaging  systems  is  usually  present  at  microwave  frequencies  only  more  so 
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because  of  the  longer  wavelength  and  the  usually  highly  thinned  nature  of 
practical  holographic  microwave  recording  aperture. 

Our  research  in  wavelength  and  polarization  diversity  was  initiated 
to  overcome  or  circumvent  those  limitations  on  microwave  imaging.  The  main 
philosophy  behind  this  approach  stems  from  consideration  of  the  "tn^osunation 
content?'  or  '''infiosuncutLon  capacity"  of  an  electromagnetic  wavefield.  A  struc¬ 
tural  measure  of  this  information  content  is  the  number  of  dcgACC6  o 6  ^Acedom 
of  the  field  [  1  ]  —  [4 ]  which  can  be  defined  as  the  number  of  independent  real 
parameters  needed  to  describe  the  field  completely.  Practical  cost  effective 
microwave  apertures  are  expected  to  be  highly  thinned  consisting  of  a  few  tens 
of  elements  or  receiving  stations  dispersed  over  the  area  of  the  aperture. 

The  amount  of  data  collected  by  such  an  aperture  intercepting  the  object 
scattered  wavefield  is  therefore  expected  to  be  quite  limited.  For  example 
with  an  aperture  consisting  of  N=*50  receivers  operating  at  a  single  frequency 
we  can  collect  50  complex  or  100  real  data  points.  Clearly  this  is  not 
sufficient  to  reconstruct  an  image  of  a  complicated  scatterer  such  as  a 
spacecraft  or  an  aircraft.  By  recording  the  field  with  the  same  array  at 

say  M=1000  distinct  frequencies  (frequency  or  wavelength  diversity)  the  number 
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of  data  points  collected  can  be  increased  by  a  factor  of  1000  to  MxN  =  5  10 
complex  points  or  10^  real  points.  By  taking  advantage  of  the  vector  nature 
of  the  electromagnetic  field  and  measuring  two  orthogonal  components  of  the 
field  each  in  amplitude  and  phase  (which  specifies  fully  the  state  of  polar¬ 
ization)  we  can  double  the  number  of  complex  data  points  to  10^,  The  net 
number  of  degrees  of  freedom  has  now  been  increased  to  2x10^.  An  impressive 
increase  by  a  factor  of  2000  over  the  single  frequency  or  monochromatic  case. 
This  is  expected  to  lead  to  significant  improvement  in  image  quality  pro¬ 
vided  we  know  how  to  process  the  degrees  of  freedom  associated  with  wavelength 


and  polarization  diversity  and  successfully  convert  them  to  spatial  image 
detail.  It  is  not  surprising  therefore  that  certain  mammals  such  as  the 
dolphin  and  the  bat  employ  frequency  diversity  in  their  echo-location  sonar 
to  probe  their  environment  with  amazing  resolution. 

A  basic  task  in  our  research  effort  therefore  is  the  analysis  of  con¬ 
cepts  and  methods  for  transforming  the  measured  data  into  useful  information 
in  the  form  of  an  appropriately  displayed  image  of  the  scattering  (or  radiating) 
object  in  which  information  contained  in  the  polarization  and  spectral  degrees 
of  freedom  of  the  scattered  wavefields  has  been  converted  into  spatial  image 
detail.  Central  to  the  practical  implementation  of  this  concept  is  the 
ability  to  take  advantage  of  recent  and  continuing  advances  in  microwave 
technology  such  as  accurate  automated  broadband  amplitude  and  phase  measure¬ 
ment,  high  speed  microwave  A/D  conversion  and  others  that  could  be  used  in 
systems  that  cleverly  substitute  recorded  spectral  and  polarization  degrees  of 
freedom  for  spatial  degrees  of  freedom  which  are  more  costly  to  record  be¬ 
cause  they  are  associated  with  the  number  of  elements  in  the  recording  aper¬ 
ture. 

Our  findings  to  date  confirm  the  validity  of  our  approach  and  are  ex¬ 
pected  to  lead  to  a  new  generation  of  imaging  radars  capable  of  high  resolu¬ 
tion  tomographic  imagery  of  distant  objects.  Detailed  discussion  of  our 
previous  findings  is  found  in  previous  publications  (see  list  of  publications). 


2.  ACCOMPLISHMENTS 

A  significant  accomplishment  during  the  period  of  this  effort  has 
been  a  mathematical  reformulation  of  wavelength  diversity  Imaging  that 
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accounts  for  the  vector  nature  of  the  electromagnetic  field  and  provides 
therefore  insight  and  suggestions  for  ways  in  which  polarization  diversity 
can  be  utilized  to  enhance  image  quality.  This  new  analy  is  provides  also 
the  mathematical  basis  for  extending  the  concepts  of  3-D  tomographic  imagery 
by  wavelength  and  polarization  diversity,  developed  so  far  for  active  imaging, 
to  passive  imaging  of  self  emitting  thermal  and  non-thermal  objects. 

An  equally  significant  achievement  during  the  past  period  was  the  fur¬ 
ther  development  of  AMTDR  (amplitude  modulated  target  derived  reference)  con¬ 
cept  and  its  verification  by  computer  simulation.  Also  during  this  period 
reconfiguration  of  our  anechoic  chamber  measurement  system  to  enable  the  ex¬ 
perimental  study  of  polarization  diversity  and  the  AMTDR  method  has  been 
completed.  Some  of  the  required  additions  to  our  instrumentation  have  been 
received  with  the  remainder  scheduled  however  for  delivery  in  the  near 
future. 

2 . 1  Vector  Formulation  of  3-D  Tomographic  Imaging  by  Wavelength  and 

Polarization  Diversity 

The  aim  of  this  section  is  to  develop  the  basic  vector  equations 
for  the  far  field  scattered  by  an  obstacle  whose  size  is  much  larger  than 
the  wavelength  of  the  incident  illumination.  The  approach  adopted  in  this 
formulation  is  based  on  reviewing  first  the  derivation  of  the  expressions 
for  the  field  radiated  by  a  bounded  three  dimensional  sinusoidal  current 
distribution.  The  results  are  then  used  to  formulate  scattering  from  a 
perfectly  conducting  obstacle  in  terms  of  surface  current  distributions  in¬ 
duced  by  the  incident  illumination.  As  compared  to  other  treatments  of 
electromagnetic  scattering,  based  on  Greenfs  Theorem  15]— [7] ,  the  treatment 
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given  here  has  the  advantage  of  avoiding  the  complexities  of  introducing 
line  current  sources  on  the  shadow  boundary  of  an  illuminated  scatterer  that 
are  needed  in  order  to  make  the  scattering  field  obey  Maxwell ! s  equations* 

We  start  from  Maxwell's  equations  for  harmonic  time  dependence  of  the 
fields  in  which  the  harmonic  time  dependence  e^^  has  been  suppressed, 


7  x  E  = 

(1) 

V  X  H  =  JuD  +  Jtotal 

(2) 

where 


^total 


J  +  J 


c 


(3) 


is  the  total  current  density  made  up  of  the  source  current  density  J  and 
the  conduction  current  density  .  The  four  field  vectors  in  eqs.  (1) 
and  (2)  and  the  conduction  current  density  are  related  by  the  constitutive 
relations, 


D  =  eE;  B  =  ]jH;  *  aE 


(4) 


where  e,  y  and  a  are  respectively  the  permittivity,  the  permeability  and 
conductivity  of  the  medium  which  is  assumed  to  be  homogeneous  and  isotropic. 
With  eqs.  (3)  and  (4)  we  can  rewrite  eq.  (2)  in  the  form. 


7  x  H  *  E  +  J 


where 


(5) 


£ 


c 


(6) 


is  the  complex  permittivity  of  the  medium. 


Solutions  for  the  electric  field  E  can  be  obtained  from  eqs.  (1)  and 
(5)  in  terms  of  the  vector  potential  A  in  the  well  known  form  [6], 

E  =  -jo)A  +  — -  V  (7- A)  (7) 

with  A  defined  by  the  relation 

B  -  V  x  A  (8) 

which  satisfies  the  third  Maxwell  equation 

V-B  =  0  •  (9) 

Defined  in  this  fashion  the  vector  potential  satisfies  the  inhomogeneous 
wave  equation, 

V2!  +  k2A  =  -uJ  (10) 

2  2 

with  k  =  cj  uec  being  the  complex  wavenumber  or  propagation  constant. 

Given  a  sinusoidal  current  source  distribution  J  bounded  in  a  volume 
V T  as  shown  in  Fig.  1,  the  solution  of  eq .  (1)  is, 

A(r)  =  fe  1  J(r')  G(r,7')  dv'  (11) 

v' 


where,  r*  and  r  are  the  position  vectors  of  the  source  point  and  field 
point  or  observation  point  respectively  and, 


G(r,r>) 


e-Jfclr-r'l 

Ir-r'l 


e 


-jkR 


R 


(12) 


is  the  appropriate  Green’s  function. 


SOURCE 

POINT 


S 


FIELD 

POINT 


Fig.  1.  Geometry  for  computing  the  electromagnetic 

field  of  a  bounded  three  dimensional  current 
distribution. 


Noting  that  the  V  operator  in  eq.  (7)  designates  differentiations  with 
respect  to  field  point  coordinates  we  can  substitute  eq.  (11)  in  eq.  (7)  to 
obtain. 


E(r)  -  /  Iff  J(r')«.,r')  + -fer 

V  c 


V  V-J(r')G(r,r')}  dV'  (13) 


consider  now  the  vector  expansion, 

V*JG  *  G  V*J  +  J*VG  O-4) 

Because  7(7')  is  not  a  function  of  observation  point  coordinates,  V*J=0. 
Therefore, 

7-JG  -  J-VG  (15> 

Consequently  in  eq.  (13), 

V(V'JG)  -  7(J*VG)  *  V*J  VG  +  J  V  VG  *  J-WG  (16) 

and  eq.  (13)  can  be  rewritten  as, 

1(7)  -  /  {  j(7*)g(7,7')  +  4^-  J(7')-vvg(7,7')}  dv-  (17) 

V 


It  can  be  shown  (see  Appendix  I)  that 


VG  =  -<jk  +  |)  G  1R  (18) 

and 

J-VVG  =|(J_*lR)lR[-k2  +3  (jk  +  j)]  -  j  (jk  +  |)  J  J>  G  (19) 

_  _  _  —  D  y  — T  *  — 

where  R  *  r-r'  and  1_  ■  —  ■  -  is  a  unit  vector  in  the  R  direction 

R  R  | - ,i 

I  r-r  I 

In  view  of  eqs.  (18)  and  (19)  we  can  cast  eq.  (17)  in  the  form, 

i(7)  *  zf*  /  1(7>)  +  [<**’>  *V  tr  [-k2  +  3  (jk  +  b 

V' 

'  R  (jk  +  R}  T(7’)]}'  dV'  (20) 

This  solution  for  the  electric  field  applies  to  the  source  free  region 
extetnaZ  to  the  volume  V*  of  Fig.  1.  It  is  customary  to  divide  this  ex¬ 
ternal  region  into  three  overlapping  zones;  namely  the  neat  zone  where  the 
value  of  R  is  such  that  no  approximation  in  the  evaluation  of  eq.  (20)  is 
permissible,  the  FteAneZ  zone  where  R  is  sufficiently  large  for  the  field 

to  be  given  to  a  good  approximation  by  terms  containing  — —  with  n  being 

Rn 

unity  or  zero  ,  and  the  6<M-  zone,  or  the  FAOunho^eA.  zone  where  R  is  such 

that  terms  in  eq.  (20)  of  the  order  of  — r—  and  higher  can  be  ignored.  In 

R 

the  Fresnel  zone  it  is  possible  to  approximate  the  Green  function  by  noting 
that  r>r'  and  therefore 
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1/2  1/2 

R  =  |r-r^  -  [(r-r ') • (r-r ') ]  -  [r  +  r'  -  2r.r'] 


-  2  1  -r' 
r 


by  the  binsmial  expansion.  Using  the  approximation  (21)  in  the  exponent 
of  Green's  function  and  R^r  in  the  denominator  we  obtain 


e-JVr  -ITT*'2  JkVT' 

G(r , r ' )  =  -2-^—  e  e 


Accordingly  in  the  Fresnel  zone  the  field  given  by  eq.  (20)  can  be  approxi¬ 
mated  by, 

-ikr  -j— r  i  kl  •  r 

E(r)  g  ~JU>4tc -  f  [J(r,)  "  (J(r,)'TR)TR]e  r  e  r  dv 


when  the  quadratic  phase  term  in  rf^  can  be  ignored  i.e.  — — -  rf^  <<  1  or, 

2r  max  * 


eq.  (23)  reduces  to  the  far  zone  field 

__  ,  -j  kr  _  _  j  kl  •  r 1 

E(r)  *  -  /  [J(r,)  "  Trl  e  r  dV’  (25) 


Where  we  have  made  use  of  the  fact  that  when  the  origin  0  is  chosen  near 
the  object  as  in  Fig.  1  or  centered  in  it  as  in  Fig.  2  the  vectors  ~r  and  H 
for  a  distant  observation  point  will  be  parallel  and  therefore  *  T  . 

The  vector  kl^  in  eq.  (25)  can  be  identified  now  as  a  scattering  vector 
kf  -  kl^  which  for  far  zone  points  is  constant  for  all  points  on  the  object. 


Fig.  2.  Scattering  geometry. 

The  integral  in  eq.  (25)  will  then  be  a  function  of  the  current  density 
distributions  J,  the  shape  of  the  object  Vf  and  the  direction  of  the 
scattered  field. 

Equation  (25)  can  now  be  used  in  the  formulation  of  the  scattering 
problem  by  referring  to  the  scattering  geometry  of  Fig.  2  and  noting  that 
for  a  perfectly  conducting  scattering  object  the  current  density  distri¬ 
bution  is  induced  by  the  incident  illumination  and  is  confined  to  the 


surface  of  the  object  in  the  form  of  a  surface  current  distribution.  The 
current  density  J(rf)  in  eq.  (25)  is  replaced  therefore  by  the  induced 
surface  current  distribution  J^r')  with  the  position  vector  F*  describing 
now  points  on  the  surface  of  the  object.  Equation  (25)  can  therefore  be 
rewritten  as, 


E(r) 


47Tr 


_  _  _  _  __  jk  .rf 

/  [J8(r')-(Js(r')-lr)  lr]  e  r 

S' 


dS' 


(26) 


When  the  object  dimensions  are  much  larger  than  the  wavelength  certain 
simplifying  assumptions  can  be  made  that  permit  expressing  the  surface 
current  density  in  terms  of  the  incident  field.  Under  these  simplifying 
assumption  that  have  come  to  be  known  collectively  as  the  pkyA'ic&t  Op£ic& 

OK  KJUuichjOii  approximation  it  is  assumed: 

(a)  The  induced  surface  currents  are  confined  to  the  illuminated  part 
of  the  surface  of  the  object.  The  effect  of  surface  currents  induced  by 

waves  creeping  into  the  shadowed  surface  of  the  object  can  be  ignored  when 
calculating  the  far  field  at  points  facing  the  illuminated  side  of  the 
object. 

(b)  When  the  radii  of  curvature  at  all  points  of  the  scatterer  are 
much  larger  than  the  wavelength  the  relation  between  the  reflected  and 
incident  fields  can  be  obtained  from  elementary  considerations  of  plane  wave 
reflection  from  a  fictitious  infinite  conducting  plane  tangent  to  the 
surface  of  the  scatterer  at  each  point.  This  is  known  as  the  "tangent  plane 
approximation"  [5]  and  under  it  we  can  express  the  induced  surface  current 
by  the  well  known  formula, 
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—  — i  ~~t 

to  the  surface  and  H  »  H  +  H  is  the  total  magnetic  field  intensity  at 

— i  — r 

the  illuminated  surface  with  H  and  H  being  the  incident  and  reflected 
magnetic  field  intensities  respectively.  The  last  equality  in  (27)  is 
obtained  by  recalling  that  for  a  perfect  conductor  the  tangential  components 
of  the  incident  and  reflected  magnetic  field  intensities  are  equal  i.e. 
n  x  *  n  x  H*.  With  these  pky^^icat  0p£icJ>  approximations  it  is  now 
possible  to  cast  eq.  (26)  in  a  form  useful  for  evaluating  the  far  field 
scattered  by  a  perfectly  conducting  object.  That  is, 


8(7  k  ) 


2~'^lcr  ,  _  __  __  jk  *r 1 

-  J  [J  (r’)  -  (J  (r’)-l  )1  ]  e  r  dS* 

4ttt  _  s  s  r  r 

bl 

1  (28) 


or  with  eq.  (27) , 


E(r,kr)  * 


-jkr  . 

-  /  [nxH1 


2ttt 


(nxH^)’l  1  ] 
r  r 


j  kr  *  r  ’ 


dSf 


(29) 


in  which  both  n  and  are  functions  of  rT  with  rf  defining  points  on  the 
surface  illuminated  surface  of  the  object. 

The  far  zone  magnetic  field  intensity  corresponding  to  this  electric 
field  is  given  by. 


H(r,kr)  -  -£■  lr  x  E(r,k.) 


Jke 


-jkr 


2iTr 


/  (n  x  H  )  x  e  r 
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jk  -r’ 


dS’ 


(30) 


These  equations  are  in  fact  identical  to  those  obtained  elsewhere  [5], [6], 
[7]  in  terms  of  surface  integrals  over  using  a  G^een'4  Theorem  fiosunu- 
ZcutLon.  The  GAeen'4  Tktosi&n  ^OAmuZcution  of  scattering  requires  however 
the  introduction  of  a  contour  integral  about  the  shadow  line  defining  the 
boundary  on  the  object  surface  separating  and  (see  Fig.  2)  to  account 
for  the  field  discontinuity  in  passing  through  the  shadow  line  in  order  to 
make  the  resulting  far  zone  fields  obey  Maxwell’s  equations.  It  is  seen 
that  in  the  treatment  given  here  this  procedure  is  not  necessary  because 
our  formulation  is  in  terms  of  bounded  induced  current  densities. 

We  introduce  next  a  vector  characteristic  object  function. 


(  J  (rf)  -  (J  (r^-l,-)  1  .  r f  on  S. 

—  — ,  v  1  s  s  r  1 

0(r  )  =  J 

J  0  . elsewhere 


(31) 


which  is  dependent  on  the  shape  of  the  surface  and  the  surface  current 
density  distribution  upon  it.  This  allows  rewriting  eq.  (28)  in  the  form. 


E(r,kr) 


-jkr 


4-nr 


/  O(r') 


-Jk_ 


dr' 


(32) 


where  we  have  replaced  the  differential  element  of  volume  dV'  in  r' -space 
by  the  equivalent  notation  dr'  and  where  the  integral. 


jk  *r' 

T(kr)  -  /  0(r’)e  r  dr'  (33) 

extends  now  over  all  space  and  is  therefore  recognized  as  the  three- 
dimensional  Fourier  transform  of  the  characteristic  function.  Accordingly 
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one  can  reconstruct  the  characteristic  function  via  the  inverse  transform 


where  r^  is  a  three  dimensional  position  vector  in  Image  space  and  dk^ 
is  the  differential  element  of  volume  in  the  Fourier  space.  It  is 
possible  therefore  in  accordance  to  this  result  to  reconstruct  the  charac¬ 
teristic  function  or  an  "image"  of  the  scatterer  by  measurement  of  E(r,k^) 
for  a  range  of  receiver  positions  r  and  wavelengths  followed  by  correcting 

the  measured  data  for  the  range  dependent  phase  by  multiplication  by  the 

4TTr  1  kx 

range  dependent  factor  e~  before  carrying  out  the  inverse  transform. 
It  is  evident  that  this  procedure  calls  for  the  use  of  polarization  dis¬ 
criminating  far  field  measurements  and  a  means  for  determining  the  object 
range  r  to  the  various  observation  points  or  receiver  locations  at  which 
E(r,k^)  is  measured. 

Of  particular  interest  in  practice  is  the  case  when  the  scattering 
object  is  illuminated  by  a  plane  wave  of  a  specific  polarization  produced 


by  a  distant  transmitter.  The  incident  magnetic  field  intensity  at  a  point 
on  the  illuminated  face  of  the  object  can  then  be  expressed  as. 


H1(r  ’)  -He  1 


where  »  kl^  is  the  propagation  vector  of  the  incident  plane  wave  and 


Here  and  are  two  orthogonal  unit  vectors  and  H*  and  H*  are  complex 
numbers.  Equations  (35)  and  (36)  can  be  used  to  represent  any  state  of 
polarization.  For  example  if  and  have  the  same  phase  the  incident 
illumination  is  linearly  polarized  with  H*  in  eq.  (35)  making  an  angle 
9  =*  tan  ^(H^/H^)  with  the  direction  of  h^  and  possessing  a  magnitude 
Hi  =»  {(H*)2  +  (H^)2}  1^2.  If  H*  and  H*  differ  in  phase  the  wave  is 
elliptically  polarized.  Circular  polarization  occurs  when  and 

their  phase  difference  is  +  tt/2.  The  tip  of  the  vector  in  any  plane 
normal  to  the  direction  of  propagation  traces  then  a  circle  with  the  sense 
of  rotation  determining  whether  the  wave  is  left  hand  or  right  hand 
circularly  polarized.  The  sense  of  rotation  or  '‘handedness"  of  the  wave 
is  determined  by  that  hand  which  when  the  thumb  is  pointed  in  the  direction 
of  propagation  has  the  fingers  curled  in  the  direction  of  rotation  of  It 
is  important  to  keep  in  mind  that  in  eq.  (35)  is  the  magnetic  field 

intensity  at  the  origin  0  of  Fig.  2  and  as  such  is  a  constant  as  far  as  the 
position  vector  variable  rf  is  concerned.  In  view  of  eqs.  (35)  and  (36)  we 
can  now  cast  eq.  (29)  in  a  form  suitable  for  plane  wave  illumination  with 
arbitrary  polarization,  namely. 


E(r,p) 


jcoy' 


rite 


2  r 


/  (n(r’)  x  H* 
V’ 


(n(r')  x  H1  •  1  }  1  ]  ejp‘r'  dV' 
o  r  r 

(37) 


where, 

p  -  k  (1^  -  lj)  (38) 

It  is  important  to  note  at  this  point  that  for  monostotic  backscattered 

field  measurements,  i.e.  when  1  *  -1.,  the  scattered  electric  field  is 

*  r  i* 

polarized  along  the  same  direction  as  the  incident  field.  Since  this  outcome 
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is  expected,  no  polarization  diversity  information  can  be  gained  about  the 
object  by  monoAtatlc  backscatter  measurements.  No  matter  what  the  state  of 
polarization  of  the  incident  field,  the  monostatic  backscatter ed  electric 
field  will  be  polarized  in  the  opposite  sense  of  the  incident  field  regardless 
of  the  shape  of  the  object.  Monostatic  backscatter  polarization  diversity 
measurements  can  not  therefore  be  employed  to  gain  information  about  the  object 
because  the  outcome  is  known.  Mi&tL&tcutLc.  measurement  of  the  reflected  far 
field  is  therefore  necessary  in  order  to  record  polarization  related  information. 
To  show  that  the  incident  and  monostatic  backscattered  fields  are  copolarized 
or  parallel  we  note  that  for  the  incident  plane  wave. 


H1  =  -  1.  x  E1 

o  n  x  o 


Therefore  in  eq.  (37) 


nxHi=^nxl,  xEi=E—  {(n#E*)  1  -  (n*l.)  E*} 

on  i  o  n  o  i  i  o 


and  therefore 


n  x  H^#l  1 
o  r  r 


{  -  (n  x  1.  x  E1) -1  }  1 

n  i  o  r  r 

(i  [U-eX-  (»-Ti)^-Tr)  Tr 


{  1  [-  (n-E*)  -  0]}  1 

n 

i  (n-E*)  T 
n  o  i 


Usings  eqs.  (40)  and  (41)  in  eq.  (37)  we  obtain 


E(r,p) 


mue-Jkr  ^ 


E1  I  (n*lJ)eJp‘r' 
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which  is  the  required  result  showing  that  when  1  ■  -  the  scattered 

field  and  the  incident  field  are  copolarized  but  in  opposite  sense.  No 
information  can  be  gained  about  the  object  by  measuring  the  state  of  polar¬ 
ization  of  the  scattered  field  because  the  integral  in  eq.  (42)  is  not  a 
function  of  the  incident  field  polarization.  Equation  (42)  forms  also  the 
basis  for  multiaspect  monostatic  wavelength  diversity  imaging  i.e,  when  the 
object  frequency  response  is  measured  with  monostatic  transmit ter /receiver  units. 

We  define  now  as  before  a  vector  characteristic  object  function. 


(  n(r')  x  Hi  -  {lT(rr)  x  Hi  •  1  }  1  . r'  on  s. 

°  o  r  r  i 

°(r’>  ”<  (43) 

I  0  .  elsewhere 


It  is  important  to  note  that  0(rf)  is  now  a  function  of  the  shape  of  the 

illuminated  surface  of  the  scatte rer  (through  the  specification  that  r*  is 

on  s^)  and  on  the  state  of  polarization  of  the  incident  illumination  though 

H1  as  described  by  eq.(36). 
o 

We  can  now  rewrite  eq.  (37)  as, 

E(r,p)  =*  ^^Tr —  ^  0(rf)  €^p’r  dr*.  (44) 

in  which  the  integral 

r(p)  -  /  0(rT)  e^p’r  dr 1  E(r,p)  (45) 

extends  over  all  r 1  space  and  is  therefore  recognized  as  the  three  dimensional 
Fourier  transform  of  the  object  characteristic  function  0(r?).  In  addition 
to  being  dependent  on  the  shape  of  the  object  the  integral  in  eq.  (45)  is 
dependent  on  the  state  of  polarization  of  the  incident  field  provided  the 
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direction  of  the  observation  points  1^  does  not  equal  -1^. 

The  object  characteristic  function  can  thus  be  determined  by  means  of 
a  three  dimensional  inverse  Fourier  transform  of  T(p)  i.e. 

--  _ 

0(r.)  -  /  f(p)  e  dp  (46-a) 

or 

E(r,p)  e  1  dp  (46-b) 

where  r^  is  a  three  dimensional  position  vector  in  M image"  space. 

Multiaspect  complex  frequency  response  measurement  with  polarization 
discrimination  can  be  used  to  effect  the  measurement  of  F(p)  over  a  finite 
region  of  p  space.  Such  measurement  carried  with  a  finite  number  of  widely 
dispersed  receivers  permit  therefore  the  accessing  of  a  finite  sampled  volume 
of  the  three  dimensional  Fourier  space  or  p  space  of  the  scattering  object. 

The  size  of  the  volume  in  p  space  in  which  the  acquired  F(p)  data  is  stored 
depends  on  the  receiver  array  size  and  geometry  relative  to  the  scatterer 
and  on  the  spectral  width  covered. 

Equations  (37)  and  (45)  are  the  main  results  of  this  analysis.  They 
form  the  basis  for  three-dimensional  imaging  of  distant  perfectly  conducting 
scatterers  by  wavelength  and  polarization  diversity.  Equation  (37)  will  also 
serve  as  the  point  of  departure  for  our  proposed  study  of  3-D  tomographic 
imaging  of  noncoherent  radiating  objects  to  be  discussed  in  the  proposed 
research  section  of  this  renewal  proposal.  The  result  in  equation  (46)  shows 
that  it  is  possible  to  access  a  finite  volume  in  the  three-dimensional  Fourier 
space  of  the  object  characteristic  function  by  complex  multistatic  polarization 
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discriminating  measurement  of  the  scattered  field  E(r,p)  as  a  function  of 

wavelength  for  as  many  scattering  directions  1^  as  practically  possible. 

This  step  is  followed  by  correction  of  the  E(r,p)  data  by  multiplication 
27Tr  1  kr  —  — 

by  the  factor  (j  eJ  )  to  obtain  F(p).  This  requires  as  pointed  out 

earlier  knowledge  of  the  range  r  between  the  scatterer  and  each  observation 
point  involved  in  the  multistatic  scattered  field  measurement.  The  complex 
(amplitude  and  phase)  nature  of  the  required  field  measurements  calls  for 
the  availability  a  common  reference  signal  at  each  receiver  in  order  to 
synchronize  the  recording  array  elements  while  storage  of  the  range  corrected  data 
T(p)  in  the  appropriate  locations  in  p  space  requires  knowledge  of  the  vector 
p  =*  kCl^  -  1^)  for  each  receiver.  Coherent  wavelength  and  polarization 
measurements  of  the  scattered  field  is  realized  in  practice  with  the  aid 
of  broadband  coherent  receivers  with  polarization  discriminating  antennas. 
Referring  to  Fig.  2  we  see  that  a  dual  polarization  antenna  at  point  P(x,y,z) 
aimed  at  the  scattering  object  can  measure  the  orthogonal  complex  field 
amplitudes  Eg  and  E^  in  the  0  and  y  directions  of  a  spherical  coordinate  sys¬ 
tem  r,9,y?  centered  at  0.  Measurement  of  Eg  and  E^  in  amplitude  and  phase 
will  completely  characterize  as  pointed  out  earlier  the  state  of  polarization 
of  the  scattered  field  and  would  double  the  amount  of  information  recorded 
as  compared  to  "scalarized"  measurements  where  the  vector  nature  of  the 
scattered  field  is  ignored  and  only  one  measurement  of  the  field  is  made. 

Once  Eg  and  E^  have  been  measured  corresponding  cartesian  components  can  be 
computed  from. 


I 
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Eg  cosy  cosa  +  sin^P 
-Eg  siny  cosa  +  cosy 
-Eg  sina 


i 


(47) 


In  this  fashion,  regardless  of  state  of  polarization,  all  field  measure¬ 
ments  irrespective  of  scattering  direction  are  reduced  each  to  three 
cartesian  components  of  E(r,p).  Accordingly  eq.  (46)  can  be  separated 
into  three  corresponding  components  O^(r^)  °z^ri^  t^e  ^ina^ 

polarization  diversity  image  being  given  by, 

2 

|0(r  .)  I  =  02  +  02  +  02  (48) 

i  x  y  y 

The  image  formulation  described  above  accounts  fully  for  the  vector 
nature  of  the  field.  Therefore  it  furnishes  a  more  precise  and  attractive 
approach  to  the  implementation  of  3-D  imaging  by  wavelength  and  polarization 
diversity  than  the  scalar  formulation  because  of  its  greater  utilization 
the  information  content  of  the  scattered  field. 

2.2  Target  Derived  Reference  Work 

The  general  aim  of  this  task  is  the  study  and  evaluation  of  data 
acquisition  methods  which  simplify  the  accessing  of  the  Fourier  space  of 
the  scatterer.  For  example  any  method  which  can  eliminate  the  need  for 
distributing  a  common  local  oscillator  to  the  various  receiver  elements  of 
the  recording  aperture  via  direct  guided  or  free  space  propagation  and 
also  eliminates  the  need  to  measure  the  range  r  from  the  object  to  each 
receiver  would  be  very  desirable  provided  that  cost  effectiveness  is  pre¬ 
served.  For  this  reason  considerable  emphasis  in  our  research  has  been 
given  to  TDR  (target  derived  reference)  techniques  that  have  the  potential 
of  providing  such  a  capability.  The  basic  philosophy  of  TDR  is  to  make 
the  target  or  scatterer  furnish  the  common  L-0  signal  to  each  receiver  to 
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serve  as  reference  for  coherent  detection  of  the  information  bearing 
scattered  field.  There  are  Important  advantages  for  TDR  such  as  dedicating 
the  entire  resolving  power  of  the  recording  aperture  to  the  target  and 
immunity  to  atmospheric  distortion  and  many  others  which  we  have  discussed 
in  previous  publications  and  reports.  The  philosophy  in  TDR  is  based  on 
the  observation  that  if  one  illuminates  the  scatterer  with  coherent  electro¬ 
magnetic  radiation  of  wavelength  sufficiently  large  compared  to  the  scatterei*s 
size  it  will  be  made  to  behave  as  a  point  scatterer  with  a  far  field  scattering 
pattern  that  is  independent  of  the  shape  of  the  object  describing  thus  nearly 
a  spherical  wavefront.  In  electromagnetic  scattering  work  this  situation  is 
referred  to  as  tow  fasizqu&ncy  or  RayZ&igk  scattering  where  the  scattered  field 
can  be  described  in  terms  of  electric  and  magnetic  dipole  moments  of  the 
scatterer  which  are  size  and  shape  dependent.  It  can  be  shown  [5]  [8]  that  the 
scattered  field  at  large  distance  r  in  the  direction  of  the  unit  vector  1^ 
when  the  scatterer  is  perfectly  conducting  and  is  very  small  compared  to  the 
wavelength  is  given  by, 

___  ,2-jkr  _ _  7  _ 

E(r,k_)  ~  —  [q-(q,1r)1r  +  ~2  n°(™  X  1^)]  (49) 

n  =4 Tjr 


where, 


i 


Here  S  is  the  total  surface  of  the  scatterer,  q  and  m  are  complex  vectors 

_ T 

describing  the  electric  and  magnetic  dipoZe.  moment  of  the  scatterer,  E 
is  the  total  electric  field  being  the  sum  of  the  incident  field  and  the 

— T 

scattered  at  the  surface  of  the  scatterer  and  H  is  the  corresponding  total 
magnetic  field  intensity.  The  second  and  third  terms  in  the  square  brackets 
in  eq.  (49)  are  complex  and  dependent  on  the  scattered  field  direction  1^. 

The  quantity  in  the  square  brackets  is  therefore  expected  to  be  a  complex 
function  dependent  on  1^  and  the  shape  of  the  object.  This  term  represents 
the  departure  of  the  scattered  far  field  from  the  perfectly  spherical  wave 
represented  by  the  first  term  of  eq.  (50).  Starting  with  eq.  (49)  it  has 
been  shown  [9]  that  the  copolarized  component  of  the  square  bracket  term  in 
eq.  (50)  is  real  for  the  monostatic  backscatter  case  1^  =  1^  regardless  of 
scatterer  shape  or  orientation.  This  result  is  intuitively  expected  to  hold 
for  small  departures  from  monostatic  backscatter.  However  as  the  bistatic 
scattering  angle  increases  a  slight  dependence  on  scatterer  shape  and 
scattering  direction  1  will  occur  resulting  in  the  slight  departures  of 
the  scattered  field  phase  from  that  of  a  perfect  spherical  wavefront.  In 
using  low  frequency  or  Rayliegh  scattering  in  a  TDR  scheme  one  would  upconvert 
the  received  low  frequency  signal  at  each  receiver  to  the  Imaging  frequency 
and  use  the  result  as  a  local  oscillator.  If  the  low  frequency  object 
illumination  is  a  phase  coherent  subharmonic  of  the  imaging  frequency  phase 
coherence  is  maintained  at  all  receiving  elements  of  the  recording  array. 

Because  the  ratio  of  the  imaging  frequency  to  the  reference  frequency 

4 

is  expected  to  be  high,  of  the  order  of  10  or  higher,  the  slightest  departure 
in  the  phase  of  the  reference  wavefront  from  perfect  spherical  will  give 
rise  to  large  phase  errors  after  upconversion  to  the  imaging  frequency  leading 
ultimately  to  unacceptable  image  distortion.  For  this  reason  and  because 
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the  low  frequency  TDR  method  would  require  in  practice  the  use  of 
separate  gear  other  than  the  imaging  frequency  gear  our  investigations  have 
concentrated  on  an  amptiXudo,  modulated  XaAgeX  dtAXvzd  method  which  we 
call  AMTDR  which  appears  to  avoid  this  problem. 

In  the  AMTDR  method  amplitude  modulated  microwave  illumination  is 
utilized.  The  modulation  frequency  co^  is  chosen  to  be  sufficiently  low 
ranging  from  a  few  tens  of  KHz  to  a  few  MHz  depending  on  the  size  of  the 
scatterer.  The  illumination  and  the  scattered  signal  received  at  each 

receiver  will  consist  therefore  of  a  microwave  carrier  at  co  with  two 

c 

closely  located  sidebands  at  co^  +  co^.  Accordingly  if  we  express  the  trans¬ 
mitted  AM  illumination  signal  in  the  usual  form, 

s^Ct)  =  A  cos  co  t  (1  +  m  cos  co  t)  (52) 

I  c  m 

where  m  is  the  modulation  index,  the  scattered  signal  received  at  a  receiver 
a  distance  r  from  the  center  of  the  scatterer  will  be 
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(53) 


where  we  have  ignored  the  phase  shift  due  to  propagation  from  the  transmitter 
to  the  center  of  the  object  that  is  common  to  all  receivers  and  where  $(oo) 
is  the  shape  and  aspect  dependent  phase  of  any  of  the  orthogonal  (0,^)  com¬ 
ponents  of  f(p)  in  eq.  (45)  and  o  is  the  K adaA  CAo 44 -4ec£o>n  of  the  scatterer. 
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We  expand  next  the  phase  4> (oj)  about  any  frequency  to^  using  the  Taylor 


ser ies. 


<J > (cjj)  =  4>(co  )  +  (to  -  0)  )  <D'(0>  )  +  -57-  (oa  -  0)  )2  <j>"  (to  )  +  . ..  (54) 


where  the  primes  designate  as  usual  differentiations  with  respect  to  oj. 
When  $(a))  is  slowly  varying  on  the  scale  of  (D^  »  w  -  0)c,  eq.  (54)  can  be 
approximated  by  the  first  two  terms  and  the  phase  of  the  two  sidebands  may 
be  approximated  by, 


4>(co  +  ca  )  *  <j)(io  )  +  a)  ’  (oj  ) 

r_  —  m  C'—  HI  C 


with  the  approximation  becoming  more  exact  as  u)  is  decreased. 

m 

By  combining  eqs.  (53)  and  (55)  we  can  write. 


s_(t)  =  Jo  Acos  [u)  (t  -  — )  +  <p(u)  )]  {  1  +  m  cos  [to  (t  -  — )+  co  (to  )  ] } 

R  ccTc  mem  o 


Demodulation  of  this  signal  by  envelope  detection  or  other  suitable 
means  yields  a  low  frequency  reference  signal 


s-  ,(t)  =  C  cos  [00  (t  -  7)  +  to  )J 

Ref  m  c  m  c 


where  C  is  a  constant.  This  signal  can  be  upconvected  in  frequency  to 
o)c  -  6  by  using  it  for  example  as  the  external  reference  in  a  microwave 
synthesizer  tuned  to  the  imaging  frequency  u j  .  In  this  manner  a  TDR  local 
oscillator  signal  for  measuring  the  phase  of  the  received  signal  of  eq.  (56) 
is  generated.  The  measured  phase  with  such  a  procedure  would  be. 
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4>  (“>  )  *  co  (t  -  tO  +  <J>(w)  -  Boo  (t  -  — )  -  6  co  <J>'  (to  ) 
me  c  c  c  m  c  m  e 


(58) 


and  since  g  -  oj  /co  , 
c  r 

<<>m(a,c)  "  *(a)c)  "  Wc  (59) 

The  measured  phase  at  any  carrier  frequency  in  the  AMTDR  scheme  is 

therefore  a  function  of  the  object  phase  and  will  therefore  yield  a  t/XSiqoA 

AsignaXuAe*.  The  measured  phase  4>m  equals  the  object  phase  <J>  minus  an  error  term 

that  is  proportional  to  the  derivative  of  the  object  phase*  We  note  Immediately 

that  at  all  points  gj  where  the  slope  <b’(w  )  =  0  <$>  —  <fc  contain  no  error 

c  c  m 

and  as  such  should  contribute  undistorted  Image  information.  Furthermore 
when  the  object  phase  is  known  or  can  be  determined  by  some  other  means  at 
one  initial  frequency  U)^  *  0)^  the  differential  eq.  (59)  can  be  solved  using 
4^  s  4*  as  initial  condition  to  yield  <j>(u))  to  be  used  in  image  retrieval. 

The  AMTDR  method  has  the  advantage  that  the  same  microwave  gear  used  to 
acquire  imaging  information  can  be  used  to  generate  the  low  frequency 
reference  of  eq.  (57)  and  that  this  signal  is  produced  by  the  same 
scattering  centers  of  the  object  that  give  rise  to  the  received  microwave 
imaging  signal. 

During  the  present  period  of  this  research  program,  we  have  completed 
the  design  of  a  unique  measurement  system  suitable  for  the  study  of  wave¬ 
length  and  polarization  diversity  imaging  employing  the  AMTDR  concept  out¬ 
lined  here.  Figures  3  and  4  show  block  diagrams  of  the  transmitter  and 
receiver  parts  of  the  system.  The  reference  frequency  in  this  system  is 
derived  from  a  high  stability  quartz  crystal  reference  oscillator  possessing 
a  short  term  stability  better  than  one  part  in  1011  for  1  sec  averaging  time. 
The  system  is  capable  of  automated  coherent  amplitude  and  phase  measurement 
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Fig.  3.  Wavelength  and  polarization  diversity  transmitter  with 
provision  for  TDR. 


COMPUTER  CONTROLLED 
POLARIZATION  SWITCHING 
FROM  MINC  LSI  11/2 


Fig.  4.  TDR  Receiver 


of  the  scattered  field  using  an  AMTDR  at  a  discrete  set  of  frequencies 

*  (N  15  +  20)  MHz  lying  in  the  (8-18)  GHz  range.  All  components  of  the 
system  that  are  not  part  of  our  present  equipment  pool  have  been  ordered. 
Some  of  these  have  already  been  received  with  the  delivery  of  the  remainder 
being  scheduled  in  the  course  of  the  next  two  months. 

Computer  simulation  studies  of  AMTDR  were  also  carried  out  during  this 

period  with  preliminary  results  indicating  that  when  the  scatterer  can  be 

considered  as  consisting  of  a  collection  of  point  scatterers  the  images 

obtained  from  the  measured  data  <b  is  identical  to  those  obtained  from  the 

m 

object  phase  <j)  provided  that  sharp  spikes  arising  from  points  where 
is  high  are  suitably  clipped.  This  is  an  important  result  the  reasons 
for  which  are  presently  under  investigation.  An  example  of  this  result 
is  given  in  Fig.  (5)  where  a  2-D  test  object  consisting  of  11  line 
scatters  distributed  in  a  format  resembling  crudely  the  outline  of  an 
aircraft  was  utilized. 

2.3  Three  Dimensional  Imaging  of  Incoherently  Radiating  Objects 

The  discussion  in  the  preceding  sections  shows  that  3-D  imaging 
of  a  smooth  perfectly  conducting  object  is  possible  based  on  multiaspect 
frequency  response  measurements  in  the  far  field  of  the  object.  If  we 
assume  the  scattering  process  to  be  linear  i.e.  the  effect  of  multiple 
reflections  between  parts  of  the  object  is  negligible  and  can  be  ignored, 
it  can  be  argued  in  accordance  to  linear  system  theory  that  the  same  infor¬ 
mation  collected  by  multiaspect  frequency  response  measurements  should  be 
obtained  by  multiaspect  impulse  response  measurements.  The  impulse  response 
and  the  frequency  response  being  Fourier  transform  pairs.  It  is  known  how¬ 
ever  that  the  impulse  response  of  a  linear  system  can  also  be  measured  with 


white  noise  excitation  and  temporal  cross-correlation  of  the  "output" 
response  with  the  "input"  excitation.  This  implies  that  it  should  be 
possible  to  use  white  noise  illumination  of  a  perfectly  conduction  scat¬ 
ter  and  cross-correlation  of  the  far  scattered  field  signal  with  the 
illumination  signal  to  access  the  3-D  Fourier  space  of  the  object  from  which 
an  image  can  be  retrieved  in  the  usual  way  by  an  inverse  Fourier  transform. 
Carrying  this  train  of  thought  one  step  further  one  can  not  avoid  raising 
the  logical  question  of  whether  or  not  it  is  possible  to  form  a  3-D  image 
of  an  incoherently  radiating  object  such  as  a  thermal  object  by  applying 
wavelength  diversity  cross-correlation  data  acquisition  and  processing 
techniques  to  its  natural  broad -band  spectrum  of  emissions.  To  answer  this 
question  we  begin  by  recalling  that  the  electrons  in  a  body  at  finite 
temperature  are  in  random  motion  with  which  we  can  associate  random  current 
distributions  that  are  therefore  spatially  uncorrelated.  These  time  varying 
current  distribution  will  give  rise  to  a  radiated  noise-like  electromagnetic 
field. 

Within  a  narrow  spectral  width  this  elec tromagnetic  field  radiated  from 
a  small  inf initisimal  segment  of  the  object  over  which  the  current  fluctuations 
are  coherent  could  be  represented  by  the  differential  form  of  eq .  (26) 
where  k  know  represents  the  mean  wave  number  in  the  narrow  spectrum  of 
radiation  considered,  and  is  the  amplitude  of  the  surface  current 
density  distribution  at  the  mean  wavenumber  k  and  is  the  propagation 
vector  in  the  direction  of  observation.  In  other  words  by  viewing  the 
incoherently  radiating  object  through  a  narrow-band  filter  we  can  regard 
it  as  a  quasimonochromatic  object  of  a  wavelength  determined  by  the 
cneter  frequency  of  the  filter  with  eq.  (26)  applicable  to  any  suffi¬ 
ciently  small  element  dS1  of  its  surface  over  which  the  current 
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distribution  is  spatially  coherent.  Accordingly  the  field  at  an  observation 
point  P(x,y,z)  due  to  a  surface  element  dS 1  of  the  incoherent  object  as 
obtained  from  eq.  (26)  and  by  referring  to  Fig.  6  can  be  expressed  by, 
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The  mutual  intensity  of  such  a  field  at  two  field  points  described  by  the 
position  vectors  r^  and  in  Fig.  6  will  be. 
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(61) 


where  *  designates  the  complex  conjugate  and  n  =  is  the  intrinsic 

impedance  of  the  medium  surrounding  the  body.  Letting  dSf  *  d£*  x  d£*  and 
noting  that  Js(r,)d£t  *  is(r*)  and  Js(rf)  dS1  =  is(rT)  d£f  where  d£f  is  a 
differential  vector  in  the  direction  of  Jg  and  ig  is  the  current  in  a 
differential  strip  of  width  d£f  we  can  rewrite  eq .  (61)  as, 
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Or  finally, 

Jl?  *  4  nk2  “ — “ -  ig (r  1 )  [1  -  cos2  9^(r ’)  -  cos2  02(rf) 


+  cos  G^(r*)  cos  02(rl)  cos0  ] 


where  the  vector 


p  =  k  -  k  =  k  (  1  -1) 

rl  r2  rl  r2 


shown  in  Fig.  6  has  a  length 


.  e 

p  =  2k  sin  — 


and  0  0  are  the  angles  between  d£*  and  1  and  1  respectively  which 

1  2  R^  R0 

are  taken  to  be  parallel  to  1  and  1  for  distant  observation  points  for 

ri  r2 

which  and  r2  are  much  larger  than  the  extent  of  the  object  as  indicated 

in  Fig.  6.  The  angle  0  is  the  angle  between  1  and  1 

rl  r2 

Because  of  the  spatially  incoherent  nature  of  the  random  surface  current 
distribution  of  the  object,  the  cross-correlation  function  or  mutual  coh&ldnce. 
fiunctton  ri2  of  the  fields  at  and  due  to  the  entire  object  ca^ 
obtained  by  integrating  eq.  (63)  over  the  surface  area  commonly  visible 
from  points  P^  and  P^.  Accordingly, 
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12  — 

Noting  that  —  ni^(r?)  has  units  of  watts/unit  area  we  define  a  cka/tcic£zAA^>£ic 

ftuncXion, 


^  \  ni2(rf)  [1  -  cos2  6, (r1)  -  cos2  0  (r1)  +  cos  0.  (r1) 
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elsewhere  (67) 


that  is  proportional  to  the  intensity  distribution  of  the  object.  Equation 
(66)  can  now  be  expressed  as, 
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(68) 


where  dr 1  designates  a  differential  element  of  volume  in  r1  -  space  and  the 
integral  extends  now  over  all  rf  and  therefore  is  recognized  as  the  3-D 
Fourier  transform  of  the  characteristic  function  I(rf)  that  is  proportional 
to  the  intensity  distribution  of  the  object  multiplied  by  an  angle  dependent 
obliquity  factor  as  defined  in  eq.  (67).  This  important  result  can  be  re¬ 
garded  as  a  3-D  generalization  of  the  Van  Citter t-Zernike  theorem  [10] .  It 
forms  the  basis  of  our  proposed  study  of  3-D  tomographic  imaging  of  in¬ 
coherent  objects  by  wavelength  diversity.  Multiaspect  spectrally  selective 
cross-correlation  measurement  of  the  far  field  of  an  incoherently  radiating 
objects  (e.g.  thermal  or  blackbody  radiators)  enables  accessing  a  finite 
volume  in  the  3-D  p-space  or  Fourier  space  of  the  choJutcXMA^ t,Lc  fauncJxon 
I(rf)  which  is  related  to  the  shape  and  intensity  or  bAsigkfneAA  distribution 
of  that  portion  of  the  object  commonly  visible  from  all  points  and  P0 
used  for  data  acquisition.  The  size  and  location  of  the  accessed  Fourier 
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space  depends  on  the  spectral  width  over  which  the  measurement  is  con¬ 
ducted  and  on  the  angular  spread  of  the  set  of  observation  point  pairs 


and  P2  selected.  As  can  be  verified  by  referring  to  Fig.  6  the  Fourier 

space  volume  accessed  in  this  fashion  is  always  a  sampled  truncated  cone 

with  the  p-space  data  in  it  falling  along  radial  lines  whose  length  according 

to  eq.  (64)  is  proportional  to  the  overall  spectral  bandwidth  of  radiation 

covered  by  the  narrowband  spectrally  selective  cross-correlation  measurement 

conducted  at  each  pair  of  observation  points  P^  and  P^*  Finally  we  see  from 

eq.  (68)  that  when  P^  and  P^  are  symmetrically  located  such  that  r^  =  =  r, 

the  range  dependent  phase  term  vanishes.  The  measured  cross-correlation  func- 

k 2 

tion  r^^p)  is  then  within  a  factor  (— )  equal  to  the  Fourier  transform  of 
the  object  characteristic  function. 

2.4  Hybrid  Data  Processing  and  3-D  Image  Display 

The  aim  of  this  task  in  our  research  effort  the  investigation 
of  methods  for  efficient  3-D  image  reconstruction  and  display  given  the 
Fourier  space  data  of  a  scattering  perfectly  conducting  object  or  of  an 
incoherently  radiating  object.  Hybrid  (opto-digital)  data  processing  tech¬ 
niques  have  received  special  attention  because  they  emerge  naturally  in 
wavelength  diversity  imaging,  and  are  capable  vf  efficient  real-time  pro¬ 
cessing  of  the  relatively  large  data  volume  involved,  and  because  they  lead 
logically  to  true  3-D  image  display  capability  which  is  highly  desirable 
if  the  full  potential  of  the  human  eye/brain  system  is  to  be  effectively 
applied  to  image  identification.  Our  approach  in  this  task  is  based  on 
Fou/UeA  domain  ptio jzcJxon  the.0A.CM6  which  we  have  discussed  in  previous 
publications  (see  list  of  publications)  and  reports  [11] .  These  theorems 


enable  the  reconstruction  of  a  3-D  image  of  the  object  in  separate  parallel 
slices  from  a  series  of  weighted  Fourier  domain  projection  holograms  one 
for  each  slice,  produced  by  digital  display  from  the  same  3-D  Fourier  space 
data  manifold.  Of  particular  interest  has  been  a  method  involving  the 
virtual  optical  Fourier  transform  (VFT)  which  allows  the  display  and  viewing 
of  a  true  3-D  image  in  the  form  of  parallel  slices  that  are  displayed 
virtually  in  rapid  succession  at  diff erent  correct  depths  in  space  in  front 
of  the  viewer.  The  details  of  our  study  of  the  VFT  and  its  application  in 
such  a  3-D  display  scheme  are  found  in  [11]  (and  in  13]  of  our  list  of 
publications) .  Although  the  use  of  the  VFT  in  3-D  display  has  many  advantages, 
it  has  the  disadvantage  that  only  one  observer  can  view  the  virtual  3-D 
image  displayed  tomographically  because  he  has  to  look  at  an  axially 
oscillating  point  source  though  a  series  of  projection  holograms,  one  for 
each  slice  of  the  object,  that  are  passed  in  front  of  his  eyes  in  rapid 
succession.  This  capability  may  be  adequate  for  certain  applications  but 
is  obviously  not  satisfactory  when  the  3-D  display  is  to  serve  several 
viewers  simultaneously.  For  this  reason  we  are  studying  at  this  stage  a 
variation  of  the  VFT  scheme  which  makes  use  instead  of  a  Real  Fourier  Trans¬ 
form  (RFT)  to  stroboscopically  project  the  reconstructed  images  of  various 
slices  in  rapid  succession  on  an  axially  oscillating  projection  surface. 
Analysis  and  design  of  an  experiment  to  verify  this  concept  has  been  recently 
completed.  The  moving  projection  screen  in  this  experiment  will  be  imple¬ 
mented  by  means  of  a  rotating  wheel  with  several  specially  shaped  teeth, 
on  its  periphery  that  furnish  the  moving  projection  surface. 

Another  central  activity  in  this  task  has  been  the  study  of  the  po¬ 
tential  capabilities  of  recyclable  SLMs  (spatial  light  modulators)  in  real- 


time  3-D  image  reconstruction  and  display*  Specifically  we  are  studying 
the  feasibility  of  optically  addressing  an  available  SLM,  and  Iteck  Corp. 

PROM  (Pockles  Readout  Optical  Modulator)  from  a  high  resolution  blue 
phosphor  CRT  on  which  we  can  display  different  computer  generated  projection 
holograms  in  rapid  succession  in  an  effort  to  reconstruct  their  corresponding 
images  at  frame  rates  of  33/ sec  or  higher .  Frame  rates  of  up  to  several 
per  second  have  been  achieved  so  far  with  reasonable  image  quality.  Higher 
frame  rates  appear  to  call  for  higher  exposure  levels  than  those  feasible 
with  the  present  high  resolution  CRT  display.  One  approach  being  pursued 
is  to  intensify  the  PROM  by  coupling  it  to  a  suitable  image  intensifier. 

A  second  generation  image  intensifier  is  being  purchased  for  our  Electro- 
Optics  Laboratory  under  seperate  funding.  This  will  provide  an  opportunity 
for  studying  the  performance  of  the  intensified  PROM  in  3-D  image  recon¬ 
struction  and  display. 

2 . 5  Design  of  Roof  Experiment 

The  object  of  this  task  is  to  generate  design  criteria  for  an 
experiment  to  image  an  aircraft  using  wavelength  diversity  imaging  tech¬ 
niques.  The  proposed  location  for  such  an  experiment  is  the  roof  antenna 
site  of  the  Moore  School’s  Graduate  Research  Center  because  of  its  proximity 
to  our  digital  storage  and  computing  facility  housed  in  the  Moore  School. 
During  this  research  period  it  was  determined  that  the  most  cost-effective 
approach  to  implementing  such  an  experiment  is  to  use  a  mono  static  back- 
scatter  data  acquisition  scheme  as  depicted  in  Fig.  7  .  A  specially  modified 
tripod  will  provide  azimuth  V*  and  elevation  0  data  with  pointing  accuracy 
of  .1  degree.  The  0,  data  will  furnish  the  polar  directions  in  the 
Fourier  p-space  for  data  storage.  Fourier  space  data  acquisition  is 
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Fig.  7.  Monostatic  Wavelength  Diversity  Imaging  Arrangement 


achieved  utilizing  the  AMTDR  scheme  with  repeated  rapid  frequency  stepping 
during  manual  optical  tracking.  In  this  fashion  the  data  acquired  during 
a  straight  line  over-flight  of  an  aircraft  passing  within  range  will  fall 
on  a  planar  surface  in  p-space.  If  the  flight  path  is  not  a  straight  line 
the  p-space  surface  will  be  curved.  It  is  hoped  that  we  will  be  able  to 
acquire  p-space  data  for  as  large  an  extent  of  the  spherical  coordinate 
y*  as  practically  possible.  Collected  p-space  data  will  then  be  projected 
on  the  horizontal  p^  -  plane  to  generate  a  projection  hologram  from  which 
we  will  attempt  to  reconstruct  a  cross-sectional  image  of  a  horizontal  slice 
through  the  aircraft.  Because  of  the  mono^taXic  mode  of  data  acquisition 
this  experiment  is  not  expected  to  yield  depolarized  field  components  from 
which  polarizatio.  diversity  information  can  be  extracted.  It  would  how¬ 
ever  be  possible  to  acquire  polarization  diversity  data  by  illuminating  the  target 
with  different  polarizations,  e.g.,  circular,  linear  vertical  and  linear 

horizontal.  Prior  to  roof  installation  we  expect  to  clear  the  experiment 
with  the  FAA  and  to  evaluate  the  proposed  system  of  Fig.  7  in  our  anechoic 
chamber  making  use  of  an  aluminum  foil  plated  toy  model  of  a  B-52  aircraft. 

This  plan  awaits  delivery  of  certain  key  components  of  the  system  of  Fig.  7 
namely  the  Quartz  reference  oscillator  and  the  absorptive  pin  modulator. 

2 . 6  Improved  Measurement  System  and  Experimental  Results 

During  this  period  two  improvements  in  the  measurement  system 
installed  in  our  anechoic  chamber  which  was  described  in  a  preceeding  report 
[11]  were  made.  A  low  noise  Avantek  AWT-18Q3911,  6-18  GHz  GaAs  FET  broad¬ 
band  amplifier  was  added  as  shown  in  Fig.  8  to  increase  the  SNR  (signal-to- 
noise-ratio) .  Enhanced  SNR  is  desirable  for  increasing  accuracy  and  improve 
image  quality.  The  gain  of  40  dB  added  by  the  broadband  amplifier  was 
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Improved  measurement  system 


selected  with  the  roof  experiment  in  mind  where  the  object  would  be  at 
a  range  of  the  order  of  a  few  kilometers.  For  the  anechoic  chamber 
measurement  an  added  gain  of  20  dB  was  found  to  be  more  than  adequate  to 
produce  drastically  improved  performance  without  exceeding  the  maximum 
input  signal  specifications  of  the  Hp  3743  Reflection/Transmission  Test  Set. 
For  this  reason  a  20  dB  attenuation  pad  was  inserted  at  the  amplifier  out¬ 
put  as  shown  in  Fig.  8.  This  measurement  system  is  semi-automatic  with 
automation  (instrument  control,  object  positioning  control,  data  acquisition 
and  storage)  being  provided  by  a  Digital  Equipment  MINC  11  Computer  System 
as  detailed  in  Fig.  9.  The  second  improvement  made  in  the  measurement  system 
was  the  refinement  in  the  azimuth-rotation  turntable  positioner  used  in 
our  experiments  to  alter  the  test  object  orientation.  The  refinement  re-r 
moved  a  nearly  unnoticeable  wobble  during  rotation.  The  effect  of  these 

improvements  is  shown  in  Figs.  10  to  13.  Figure  10  is  a  comparison  of  the 
improved  system  imaging  performance  with  theoretical  performance  determined 
by  simulation  as  described  in  an  earlier  report  [11].  Excellent  agreement 
between  the  theoretical  and  experimental  projection  holograms  and  tomographic 
images  of  the  two-cylinder  test  object  used  earlier  [11]  is  seen  to  result. 

The  image  represents  a  horizontal  slice  l  v  ough  the  2-vertical  cylinders. 

This  improved  performance  is  further  detailed  in  Figs.  11  to  13.  Shown  in 
Fig.  11  is  the  (6-18)  GHz  fixed  aspect  frequency  response  of  the  two-cylinder 
test  object  described  in  [11].  The  (6-18)  GHz  frequency  range  is  presented 
in  three  frequency  bands.  For  each  band  256  complex  (amplitude  and  phase) 
data  points  were  taken  using  the  measurement  procedure  described  in  [11]. 

The  frequency  response  shown  for  vertical  linear  polarization  is  not  corrected 
for  system  response  or  room  clutter  which  are  given  in  Fig.  12.  Both  system 
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(c)  (d) 


Fig.  10.  Comparison  of  improved  imaging  system  performance 
with  theory.  Top  -  projection  holograms,  (a) 
measured  (b)  computed.  Bottom  -  optically  recon¬ 
structed  images  with  zero  order  light  removed  by 
spatial  filtering,  (c)  image  retrieved  from  (a), 
(d)  image  retrieved  from  (b) . 

4  3 


oooo*«-oi  Y.  looot  *-oi i  ioooc*ci 
:  WIT  MIA  f  HE  NM*  : 

’  6T  t :  1 1 .  M  r 

sta*t  uffl.tM;.  ia  stop  n*o.««.  t*  •  m:a  poihts  m 


01 

.a.  ^  _  XVs 

s  ;  \  /  \ 

U  f- 

i 

•L- . 

i  . .  i .  i 

T  *,  *  0<X»  •  o  1 

\  *ooo*  *01  Ti  .ooo*  *ri  i  t  ooo* 'Oi 

J  y  —  .  ... - - —  r 

l 

L_ 

l  «ooo»«oi 

t  i  i ! 

Vvkkx.oT  i.fcCKKx-ot  t.;oo«'6i 

I  »  DATA  P0IHIS 


Fig.  11.  Uncorrected  frequency  response  of  two-cvlinder  test 

object.  (a)  6-10  GHz,  (b)  10-14  GHz,  (c)  14-18  GHz, 
Polarization:  Linear  vertical.  Amplitude  in  dB,  phase 
in  degrees. 
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Fig.  12.  Lef t: uncorrected  clutter  data  for  anechoie  chamber. 

Subtract  60  dB  from  readings  to  obtain  clutter  level. 
Right? corrected  reference  target  signal  taken  as  system 
response.  (a)  6-10  GHz,  (b)  10-14  GHz,  (c)  14-18  GHz, 
Polarization:  Linear  vertical,  amplitude  in  dB,  phase 
in  degrees. 
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Fig,  13.  Frequency  response  of  two  cylinder  test  object;  Left  - 

measured  data  of  Fig.  11  corrected  for  clutter  and  system 
response.  Right  -  computed  frequency  response.  (a)  6-10 
GHz,  (b)  10-14  GHz,  (c)  14-18  GHz,  Polarization;  Linear 
vertical,  amplitude  in  dB,  phase  in  degrees. 


response  and  clutter  given  in  Fig.  12  were  obtained  using  the  procedures 
described  before  [11]  .  Removal  of  system  response  and  clutter  errors 
from  the  uncorrected  object  response  yields  the  corrected  response  shown 
in  the  left  part  of  Fig.  13.  Comparison  of  this  corrected  response  with  the 
computed  frequency  response  shown  in  the  right  of  Fig.  13  indicates  very 

TT 

good  agreement.  The  constant  phase  difference  of  about  —  seen  is  immaterial 
and  is  attributed  to  a  slight  error  in  the  distances  between  the  phase 
centers  of  the  transmitting  and  receiving  antennas  and  the  object  assumed 
in  the  computation,  (see  ref.  11  for  details  of  computation). 

3.  CONCLUSIONS 

We  have  presented  a  physical  optics,  vector  formulation  of  three  dimen¬ 
sional  tomographic  imaging  from  far  field  scattered  data  involving  wave¬ 
length  and  polarization  diversity.  To  avoid  unnecessary  complications, 
the  treatment  specifically  addresses  microwave  and  millimeter  wave  imaging 
of  perfectly  conducting  objects  which  are  characterized  in  the  analysis  by 
surface  current  distributions  induced  by  incident  plane  wave  illumination. 
Such  conducting  objects  are  of  practical  interest  in  imaging  radar.  The 
analysis  could  be  extended  however  to  weakly  scattering  dielectric  objects 
in  terms  of  volume  currents  caused  by  induced  polarization  and  to  the 
imaging  of  three  dimensional  distributions  of  refractive  index  in  the 
optical  regime.  The  resulting  equations  represent  a  vector  generalization 
of  the  inverse  scattering  imaging  problem  and  multifrequency  holography. 
Their  vector  nature  permits  accounting  for  the  effect  of  field  polarization 
on  the  retrieved  image.  The  results  show  that  it  is  possible  to  doubly 
and  nonfi&dundcivitly  access  a  finite  volume  in  the  3-D  Fourier  space  of  a 


47 


scattering  object  by  one  set  of  multistatic,  far  field,  polarization 
sensitive  complex  frequency  response  measurements  that  can  be  carried  out 
by  a  small  number  of  widely  dispersed  receiving  stations  or  observation 
points  that  define  the  recording  aperture  size.  Fourier  domain  data 
collected  in  this  manner  permit,  through  the  use  of  Fourier  domain  pro¬ 
jection  theorems,  the  retrieval  of  three  dimensional  image  detail  tomo- 
graphically  either  in  parallel  slices  that  lie  normal  to  any  desired  view¬ 
ing  directions  or  in  central  or  meridional  slices  through  the  object. 

Such  3-D  tomographic  image  reconstruction  can  best  be  carried  out  in  a 
hybrid  (opto-digital)  computing  scheme  which  is  able  to  automatically 
incorporate  polarization  information  in  the  retrieved  image  doubling 
thereby  its  information  content.  The  theory  developed  shows  that  it  is 
possible  to  obtain  high  resolution  imaging  in  a  cost-effective  manner. 

This  is  achieved  by  replacing  expensive  spatial  (or  angular)  degrees  of 
freedom  that  are  determined  by  the  number  of  available  receiving  stations 
with  less  costly  spectral  and  polarization  degrees  of  freedom.  This  trade¬ 
off  feature  has  important  implications  in  microwave  imaging  of  distant 
scatterers  specially  airborne  and  space-borne  objects  since  it  can  reduce 
the  cost  of  systems  to  a  realistic  level.  The  realization  of  high 
resolution  3-D  tomographic  imaging  radars  with  such  features  requires 
however  the  capability  of  furnishing  a  common  reference  signal  to  each 
receiving  station  within  the  recording  aperture  to  enable  coherent  measure¬ 
ment  of  the  scattered  field.  For  widely  dispersed  receiver  arrays,  of  the 
type  needed  in  high  resolution  imaging  of  distant  objects,  the  distribution 
of  such  a  common  phase  reference  or  local-oscillator  signal  by  guided  (trans¬ 
mission  line)  or  by  broadcasting  a  signal  is  unrealistic  because  of 


high  anticipated  cost  and  proven  unreliability.  For  this  reason  techniques 
that  seek  to  furnish  such  a  reference  in  a  more  cost-effective  and  reliable 
fashion  are  of  primary  importance.  Among  these  are  techniques  that  logically 
seek  to  derive  the  reference  from  the  scattering  object  itself  by  low  fre¬ 
quency  illumination  where  the  wavelength  is  sufficiently  large  compared  to 
the  object  size  to  make  it  behave  as  a  point  scatterer.  A  promising  variation 
of  this  basic  approach  employs  low  frequency  amplitude  modulation  of  the 
illuminating  microwave  signal  to  indirectly  produce  a  low  frequency  reference 
illumination.  This  AMTDR  scheme  appears  to  be  able  to  furnish,  at  the  least, 
a  3-D  tasiQ&t  &'ignatuA.&  and  can  furnish  a  3-D  geometrical  image  if  the  phase 
of  the  object  ai  the  initial  frequency  can  be  determined  for  the  various 
receivers.  In  addition  to  eliminating  the  need  for  a  local-oscillator  dis¬ 
tribution  network  the  AMTDR  scheme  has  several  attractive  features  the  most 
important  of  which  are:  (a)  It  makes  the  imaging  process  immune  to  atmospheric 
distortion  because  the  reference  signal  and  the  imaging  signals  traverse  the 
same  path.  (b)  Because  the  reference  signal  is  derived  from  the  same 
scattering  centers  on  the  target  that  produce  the  imaging  signal,  the  Doppler 
frequency  shift  in  the  upconverted  reference  signal  due  to  target  motion  and 
in  the  imaging  signal  will  be  equal  and  will  therefore  cancel.  (c)  It  is 
well  known  from  considerations  of  resolution  in  Lensless  Fourier  transform 
holography  that  the  placement  of  the  reference  source  on  or  in  the  vicinity 
of  the  object  permits  dedicating  the  entire  resolving  power  of  the  recording 
aperture  to  the  object  alone  resulting  in  higher  resolution  in  a  manner 
similar  to  zoom  photography.  (d)  Because  the  frequency  components  of  the 
AMTDR  signals  lie  in  the  microwave  range,  operation  through  the  ionosphere 
is  possible. 


The  principles  of  A  and  polarization  diversity  imaging  presented 
here  are  also  applicable  to  smaller  size  microwave  and  especially  mmW 
apertures  such  as  those  that  can  be  envisioned  for  overcoming  the  short¬ 
comings  of  electro-optical  systems  in  smoke  and  dust  filled  environments. 

In  such  cases  the  receiver  arrays  would  most  probably  be  sparsly  mounted 
on  the  surface  of  a  vehicle  in  a  conformal  manner  rendering  it  thus 

.  As  long  as  the  position  vector  of  each  array  element  relative 
to  a  vehicle  mounted  illuminator  is  known  together  with  the  direction  of 
illumination  it  should  be  possible  to  specify  the  p-vector  of  each  receiver 
element  and  therefore  the  p-space  locations  for  data  storage  and  subsequent 
image  retreival.  The  size  of  vehicle  mounted  conformal  arrays 

allows  now  direct  distribution  of  a  local  reference  to  the  various  receiver 
elements  locally  via  transmission  lines  instead  of  a  TDR  system.  Conformal 
arrays  of  this  type  would  particularly  be  attractive  for  Jjnq  the 

skin  of  an  aircraft  to  furnish  it  with  a  nearly  omni-directional  seeing 
capability.  Conformal  imaging  arrays  of  much  larger  size  can  also  be  en¬ 
visioned  in  space  application  and  in  naval  applications. 

An  advantage  of  formulating  the  electromagnetic  scattering  problem 
in  terms  of  induced  surface  currents  is  the  ability  to  extend  it  to  the 
3-D  imaging  of  incoherently  emitting  objects  by  relating  incoherent  emission 
to  random  current  distributions.  This  extension  reveals  that  it  should  be 
possible  to  access  the  3-D  Fourier  space  of  the  3-D  intensity  or  brightness 
distribution  of  such  an  object  by  wavelength  selective,  cross-correlation 
measurement  of  the  far  field  of  the  object.  The  application  of  the  3-D 
hybrid  image  reconstruction  techniques  to  this  case  indicates  the  feasibility 
of  3-D  radiometry. 
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7. 


APPENDIX  I 


i 


VECTOR  EXPANSIONS  OF  VG  and  J-WG 


It  is  convenient  to  work  in  cartesian  coordinates*  The  final 
however  are  general.  We  begin  with 


VG  =(-—•  1  +  ~  1  +  X  )  — — — 

9x  x  9y  y  9z  z  R 
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and  since, 
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we  have, 
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which  is  the  required  expansion  of  VG.  It  is  an  easy  matter  to  show  and 

3R  1  3R 

worthwhile  to  note  that  since  — f  ■  -  t-(x-x')  -  -  —  with  similar 

yX  R  c)X 

*  .  c  3R  j  3R 

relations  for  — ,  and  — , 

3y  3zT 


V1  G  =  -  VG 


with  the  prime  on  the  "del”  operator  designating  differentiations  with 
respect  to  the  primed  source  coordinates. 

We  turn  next  to  J*VVG  which  we  write  using  eq.  (6)  as. 


J-VVG 
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Using  eqs.  (4)  in  eq.  (9)  we  obtain. 
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which  can  be  reduced  finally  to  the  required  expansion, 
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